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Abstract 

Resonant inelastic (Raman) light scattering off neutral GaAs quantum dots 
which contain a mean number, = 42, of electron-hole pairs is computed. 
We find Raman amplitudes corresponding to strongly collective final states 
(charge-density excitations) of similar magnitude as the amplitudes related to 
weakly collective or single-particle excitations. As a function of the incident 
laser frequency or the magnetic field, they are rapidly varying amplitudes. 
It is argued that strong Raman peaks should come out in the spin-density 
channels, not related to valence-band mixing effects in the intermediate states. 

PACS: 78.30.-j, 78.67.Hc 
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I. INTRODUCTION 

The experimental study of the inelastic (Raman) scattering of light in arrays of quantum 
dots (qdots) began a few years ago0i. These studies were aimed at investigating multipole 
excitations or spin excited states in the quantum dot, which do not leave traces in lumines- 
cence or absorption experiments. In both papers, quasi-bidimensional dots with radii around 
100 nm and nominal electron density about 8 x 10^^ cm~^ were studied. A rich spectrum 
of single-particle (SPE), and collective charge- density (CDE) and spin-density excitations 
(SDK) was observed. 

The single-particle spectra were interpreted in terms of the single-particle density of 
states in the dot0, computed in the Hartree-Fock (HF) approximation. The energy position 
and Raman strengths of CDE states in dots with 12 electrons, computed within time- 
dependent HF theory!, were shown to qualitatively agree with the experimental results. 
More consistent calculations in the framework of time-dependent density-functional theory 
(DFT) were carried out in Ref. ^, where the charge and spin dynamic structure functions 
were computed for a system of 200 electrons. The multipolarities of the observed CDE and 
SDE peaks, and the relative peak intensities, as functions of the transferred wave vector of 
the light, were reproduced. 

Let us stress that, in electron qdots, Raman processes in which the final states are SDE 
require the account for valence-band mixing in the intermediate hole state. For this reason, 
CDE peak intensities reported in Ref. ^ can not be properly compared with intensities in 
the SDE channels. The inclusion of valence-band mixing effects in the intermediate states 
of Raman processes i. to be published etewherj. 

In the present paper, we compute Raman amplitudes for neutral quantum dots, where 
the number of holes, N, equals the number of electrons in the dot. To the best of our 
knowledge, there are no similar calculations in the literature. 

Although the lifetime of the excitons may pose certain difficulties to the Raman mea- 
surements, we believe that it does not represent a real challenge to present experiments. 
With relative independence on the dot parameters, the multi-excitonic system reaches typi- 
cal densities around 5 x 10^^ pairs/cm^, a value which may be obtained by pumping a neutral 
qdot with a stationary mean-power laser. In fact, densities well above 10^^ cm~^ have been 
achieved already for a few year^l with pulsed high-power lasers. 

Due to the lack of experimental results, we will focus on the qualitative aspects following 
from our calculations. A simplified two-band model of a disk-shaped qdot with parabolic 
confinementlll is to be used. 

We shall, first, make some remarks concerning the computational method. The Random 
Phase Approximation (RPA) is the common frame, as in the cited papers. We adopt the 
wave-function approach of Nuclear Physics!, and construct RPA approximations to the 
wave functions of both final and intermediate states. Coulomb interactions and collective 
effects are exactly accounted for within the RPA, even for the intermediate states entering 
the Raman amplitudes, which are states with + 1 e-h pairs. Corrections to the RPA 
functions, such as the e-h pairing; correlations, could, in principle, be included by means 
of the quasiparticle RPA schemeS. The electron-radiation (e-r) interaction hamiltonian is 
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written in second quantization in the basis of HF single-particle states. Our treatment of 
Raman scattering follows the lines of Refs. ^ and 10, in the sense that the e-r interaction 
causes transitions between multi-excitonic states. Thus, energy denominators contain A^- 
pair instead of single-particle energies. 

Concerning the numerical results, there are a few points to stress. First, the absorption 
threshold, or the frequency for which extreme resonance is achieved in Raman scattering, 
grows at a rate of 0.3 meV per pair added to the dot. This is an indirect way of determining 
the mean number of pairs in the dot. Second, Raman peaks in quadrupole channels are 
1/10 of monopole peaks at momentum transfer around 0.8 x 10^ cm~^. Next, intensities 
corresponding to weakly collective or SPE are comparable in magnitude to the strongest 
CDE peaks, and vary very rapidly with the magnetic field or the frequency of the incident 
laser. Thus, our calculated spectra resemble more the complex spectra of quantum wells in 
strong magnetic fieldjni, and differ from the smooth experimental curves obtained in Ref. |^ 
for the pure electronic strong-confinement qdot. 

A last point which deserves attention is the fact that Raman scattering in SDK channels 
do not require mixing of hole bands. Thus SDK peaks should be observed in any polarization 
of the scattered light. 

The present calculation of Raman cross sections completes a series of papers on the 
optical properties of the A^-pair system in a qdot. The dominance of a giant dipole resonance 
in the infrared absorption, which position scales as A^^/^, was shown in Ref. [l^. This 
resonance could be studied through the modulations of the photoluminescence caused by 
an infrared source, in the same wav as the infrared excitations in the 'W-electron plus one 
hole" system are presently studiedli^. On the other hand, the position and intensity of the 



coherent magneto-luminescence peak were computed in Ref. 14 



The plan of the paper is as follows. The basic expressions for Raman cross sections along 
with HF, RPA and particle-particle Tamm-Dancoff approximations (pp-TDA) are presented 
in Section II. The formalism is well established in the Nuclear Physics contexti. We underline 
in that section only the main points for the sake of completeness. Computed ground-state 
properties, the multipole excitations and their strength variations with magnetic field, the 
band gap renormalization as a function of the number of pairs, and the Raman cross sections 
are given in Section III. Final remarks are presented at the end. 



II. THE BASICS 

The resonant inelastic (Raman) light scattering off a neutral quantum dot containing a 
mean number, A^, of electron-hole pairs is schematically represented in Fig. The energy 
of the incident photon, hui, is supposed to be close to the band gap energy, Egap. In Fig. |l| 
(a), the final state has the same spin quantum numbers as the initial (ground) state of the 
A/^-pair system. It is, in general, a charge- density excitation (CDE). The electron and hole 
spins are represented by arrows. Incident and emitted photons are drawn as wavy lines. 
Additionally, there are also Raman processes in which the final states involve changes in 
the spin quantum numbers. These states will be called spin-density excitations (SDE). We 
have represented in Fig. |I] (b) a situation in which the total electron and total hole spin 
projections experience changes AS^e = —1 and ^S^h = +1, respectively. 

The amplitude for the Raman processes depicted in Fig. |I] (a) is given by: 
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^^^^ _^ ^Nf + l ^{f\H+_^\tnt)N+, M^,{int\H-_^\i)N^. 



int 



huji - {Eint - Ei) + iV 



where A^j, Nf are the mean number of photons in the initial and final states. For spon- 
taneous Raman scattering, Nf = 0. The sum runs over intermediate states with + 1 
pairs and the appropriate quantum numbers. F is the lifetime broadening. We will take 
it phenomenologically as F = 0.5 meV.i The resonance condition means that the leading 
contribution to (0) comes from intermediate states satisfying 

huJi ^ Eint - Ei, (2) 

i. e. the incident photon has nearly the same energy as the jump in energy from the initial 
to the intermediate quantum-dot state. 

H~_r is the interaction hamiltonian corresponding to the annihilation of a photon and 
creation of a new electron-hole pair. Its matrix element is written as 



N+l 



y,N+i{int\e^Ji\\i)N. (3) 



In the normalization factor entering Eq. (^), e is the electron charge, mo is the electron 
mass in vacuum, V is the sample volume, and r]i is the refraction index at frequency cUj. On 
the other hand, the first factor in the sum comes from the (band) spin quantum numbers 
of the initial and final states in the interband transition. We give its detailed expression in 
Appendix ^. Si is the polarization vector of the incident light. 

The next factor in the sum depends on the orbital (envelope) one-particle wave functions. 
As a basis for the one-particle states, we will use the Hartree-Fock (HF) set, ipae and (p^h 
for electrons and holes respectively, is the wave vector of the incident light. Detailed 
expressions for the orbital factor are given also in Appendix 0. 

Finally, we will compute the last factor in (|^) by means of the so called particle-particle 
Tamm-Dancoff Approximation (pp-TDA) formalism!, to be briefiy described below. 

The second matrix element entering Eq. (|l]) is written as 



N 



{f\Ht.Aint)N+i 



2'Kh 

7 



mo\ Vuffi] i 



E (^7 - P.;) (/ e-^'^^ >.e(rl0,.(rldV) 

XAr(/|/l^e„|mt)Ar+i. (4) 



Its computation involves similar band and orbital factors (see Appendix 0). The last factor 
requires - besides the intermediate (A^ + l)-pair states, to be obtained from the pp-TDA 
equations - the knowledge of the A^-pair excited states, which will be obtained from the 
ordinary or "particle-hole" RPA formalism!. The starting point for both RPA and pp-TDA 
schemes are the HF single-particle states. 
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A. The HF equations 



We take for the HF single-particle functions the following ansatz: 

0f^ = \/|sin(^)EC^i^^X.(ni), (5) 

where < z < L, L = 12 nm is the height of our disk-shaped qdot, fji is the projection of 
r onto the xy plane, and Xs are the two-dimensional (2D) oscillator wave functions, given 
elsewherelli. The Ca,s coefficients are obtained from the equations 

E 



12U15. 



j:{^,u\i/r\t,v)c!;^^c!;^^ = (6) 

and a similar set of equations for the C^^. Notice that the E^^^ are 2D oscillator energies: 

= n^Lui + Lu^jA {2h + |/.| + 1} + ^Is + gefisBS:, (7) 
4? = — nsju;i + u^j4 {2K + |/.| + 1} - ^^Is - 9h^^bBS^, (8) 



where hujQ = 3 meV is the in-plane confinement, Uc = eB/{mec) is the electron cyclotronic 
frequency, and ge, gn are Lande factors. We took parameters appropriate for GaAs: nie = 
0.067mo, i. e. hUc/B = 1.728 meV/Teslas, the ratio of in-plane masses is me/mh = 
0.067/0.11, and gefJ^B = —0.0173 meV/Teslas, ghfJ^s = 0.0296 meV/Teslas. In the gf-factors, 
the effect of qdot height was approximately accounted foJii. 

On the other hand, the {s,u\l/r\t,v) are Coulomb matrix elements taken over 2D oscil- 
lator wave functions0, and the strength /? is given by 0.8 e^/{Klo), where: 



In 



h , , 

(9) 



is the unit of length, k = 12.5 is the dielectric constant, and the 0.8 coefficient takes care 
approximately of the effect on Coulomb interaction of awaging over the z coordinatdi. 

Equations (|]) are solved iteratively. We start by occupying the lowest oscillator shells, 
construct the matrix inside brackets in (P), and iterate until convergence is reached. The 
occupation of HF levels is actualized after every 10 steps in accordance to the current values 
of the HF energies. 15 oscillator shells are used in the calculations, i. e. a total of 240 2D 
oscillator states. 
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B. The RPA equations 

In the RPaI, we allow a general correlated ground state, \RPA), and the excited states 
are looked for in the form 

^ = Q^RPA), (10) 
where the operator for CDE states is given by the expression 

QhoE = T^iX^xeiex + X^.hlh, - Y^A^. - Yth\K). (11) 

The index A runs over occupied HF states, and a runs over unoccupied states. The X, Y 
coefficients are nonzero only for transitions respecting the selection rules, i.e. the spin pro- 
jection doesn't change, and the change in angular momentum is fixed (a given multipolarity 
of the excitation). These coefficients satisfy the equations: 

T,fl 

E {Bt,r,X% + Bt,^Xr, + ^t,,rY',r + ^a'^.^^. } = "^^^t, (12) 



in which fiVt is the excitation energy, r and /i are indexes similar to a and A, respectively, 
and the A and B matrices are given byEl'Ill: 

Al\^^ = {E,, - E,x)6„Jxf. + P {{(T, fi\l/r\X, r) - {a, /i|l/r|r. A)) , 

A:l^^ = -/3{a,f^\l/r\X,T), 

BTx,,r = /3((a,r|l/r|A,/x) - (a, r|l/r|/i. A)), 

BX, = -P{a,T\l/r\\,fi). (13) 

Notice, for example, that in Al\ .^^, a and A are electronic HF states, and t, fi - hole states. 
A^'^ has formally the same expression as A'^'^, A^'^ the same as A'^^, etc. Let us stress also 



that Coulomb matrix elements over HF states enter the RPA equations ([121) , they can be 
computed from the matrix elements over oscillator states by means of the expansions (^. 
Usually, positive (physical) and negative (unphysical) excitation energies come from (|12|). 
The physical solutions annihilate the RPA ground state 

Q\RPA) = 0, (14) 

and satisfy the normalization condition 

1 = Eii^^Ar + i^.\r - \Y^' - \Yt\'}. (15) 
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To evaluate the collective character of a state we compute the matrix elements of the 
multipole operators: {"^iDilRPA). Collective states give significant transition strengths, 
whereas single-particle excitations give practically zero matrix elements. The multipole 
operator, Di, is defined as: 



(16) 



a,7 



where 



a7 



(«|p2|7); / = 0. 



(17) 



p and 9 are polar coordinates in the xy plane. Detailed expressions for c?'^^ are given in 
Appendix ^. 

Multipole matrix elements are computed from the RPA amplitudes in the following way: 



They fulfill the energy-weighted sum rulesi 



(18) 



1 



+ — E (AI(r^)'^'-^IA) 



(19) 



for / 7^ 0, and 



Y^hQ^ mDo\RPA)\' = 2h'e'{ Y: (A| — |A)+ Yl (A| — |A) 



r 



(20) 



for / = 0. The /i^'s are Fermi levels. Thus, X < fip means that the sum runs over occupied 
HF states. Explicit evaluation of the r. h. s. of ( P^pOD is done in Appendix 

Spin excitations can also be built on within the RPA formalism. For example, a state 
with ASez = 1, AShz = can be obtained from a like (|n|) with only electron operators, 
such that the transitions satisfy the spin selection rule. It should be noticed, however, that 
the simple combination of one-particle excitations in Eq. (0) does not allow us to construct 
"2j9 — 2h" excited states with ASez = 1, AShz = — 1, for example, entering the final states 
of Raman SDE processes. 



C. pp-TDA 

The pp-TDA scheme allows us to build up states with 2A^ + 2 particles starting from 
the the ground state of the A^-pair system!. Notice that there are 12 possibilities for the 
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added pair of particles. We can add, for example, an e-e pair with various spin orientations. 
In this subsection, we focus on the situations where an optically created e-h pair is added. 
That is, only the following two possibilities are considered: ethj or e|h|. 
The Q'^ operator, analogous to Eq. (pH]) , is written in the following form: 

gt = ^V;(f+^)et/,t, (21) 

(T,r 

where, as before, a and r label states above the Fermi levels. acts on the RPA ground 
state to produce states with + 1 pairs. The V coefficients satisfy the equations: 

m -K- EMr^'^ = E AVr\a\ r')ri,7^\ (22) 

<7',t' 

The quantity hVL gives the excitation energy, measured with respect to the \RPA) A^-pair 
state: E{N + 1) - Erpa{N). 

RPA and pp-TDA excitations energies and coefficients X, and V are to be used in the 
computation of Raman amplitudes. 



D. Raman scattering in CDE channels 

The inelastic scattering of light, schematically represented in Fig. |I] (a), is characterized 
by a Raman shift: huji — hujf = Ef — Ei = hQf. The amplitude for the process is given by 
Eq. (|I|). This amplitude will depend on the scattering angles. 

To state a convention, the dot plane will define the xy plane, and the magnetic field - 
the positive 2;-axis. The incident light comes from the z < subspace, forming an angle (pi 
with the z axis. The g^, B pair of vectors define the xz plane, i. e. the projection qi\\ is 
oriented along the positive x axis. The emitted light goes back to the z < subspace. It 
is characterized by angles 0/ with the z-axis, and 6'/ in the xy plane. We will take, for the 
incident light: 

gill = qi sin 

s,y = l, (23) 

whereas for the scattered light: g^y = g^sin^j. We will distinguish two situations: (i) The 
"parallel" light polarization, in which: 

efcc = -sm9f, 

efy = cose f, (24) 
and (ii) The "perpendicular" light polarization, where: 

Efx = COS (f) f COS 9 f , 

e fy = COS (pf sinOf. (25) 
Below, we give an explicit expression for the matrix element of the H~_^ hamiltonian: 
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MM^nt\H;_^\t)N = —J^^T^banS^l x ortoa/W x V^^+'> . (26) 
mo V VuJi7]f ^ ' ' ' 

The band and orbital factors are evaluated in Appendix 0. The energy denominator in 
the scattering amplitude will be written in the form: 

huji - {Eint - Ei) = huji - Egap - hn\^^^^\ (27) 

where hflint is the eigenvalue coming from the pp-TDA equations, and Egap = 1560 meV is 
a nominal band gap. 

The matrix element of the H^_j. operator, has a bit more cumbersome expression: 



N{fm_,\znt) 



N+l 



mo \ VufT]} ^ 



^ A ^2 Y.Y.K^y) bandit X orUtali^^ x Vj^^'^X^r,. (28) 

The interpretation is, however, straightforward. Let us take the first term. The first sum 
runs over e-h states {(t,t), both above the Fermi levels, entering the pp-TDA function 
\int) jv+i. V^^^"^^-* are the corresponding coefficients. The second sum represents the electronic 
excitation part of the RPA function |/)iv. A is an electronic state below ^p, and a' an state 
above fip. The transition from \int)N+i to \ f)N is caused by a pair of annihilation operators 
he. It is evident that the subindexes should be hrCx, and thus cr' = a. 

Amplitudes for backscattering processes, in which (j)f = (pi, and 9f = n, will be computed. 



III. RESULTS 
A. Properties of the HF ground state 

We give in this sub-section a few qualitative results that follow from the HF calculations. 

We show in Fig. ^ a subset of the HF single-particle levels at i? = 1 T. Fermi energies are 
represented as dotted lines. Apart from an overall downward shift, we observe only a slight 
deformation of free oscillator shells due to Coulomb interactions. The Zeeman splitting is not 
resolved in the figure scale, thus spin-up and -down levels are simultaneously occupied. As 
a result, total electron and hole spins remain equal to zero when the magnetic field is varied 
between zero and 2 T. In fact, very low spin polarizations persist up to higher magnetic 
fields, of the order of 20 T.0 The total ground-state angular momentum is also zero in this 
magnetic field range, and persists up to very high B, as a prelude to the formation of e-h 
pairs in zero relative angular momentum states, which maximize Coulomb attraction. 

Let us stress also that Fig. ^ qualitatively predicts that single-particle excitations (SPE) 
with ALz = ±2 ("quadrupole"-like, represented by arrows in the figure) are lower in energy 
than "monopole" (AL^ = 0) or "dipole" (AL^ = ±1) excitations at i? = 1 T. This fact is 
corroborated by the RPA calculations, see below. The energetic cost of adding an e-h pair 
with le + Ih = is, according to Fig. ||, around 15 meV (plus Egap). A value confirmed by 
the TDA results. 
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Electron and hole densities at B = and 2 T are drawn in Fig. ^ The small differences 
between both densities are due to the differences between electron and hole in-plane masses. 
The maximum value, around 7 x 10^^ pairs/cm^, is typical of excitonic systems. Notice 
also that density oscillations, related to shell-filling effects, are smoothed as B is increased. 
This fact is due to the increasing occupations of states in the first Landau level, which wave 
functions have no radial nodes. 

B. Multipole excitations and renormalization of the absorption edge 

CDE of various multipolarities in the A^-pair system are obtained from the RPA com- 
putations. We show in Fig. ^ the monopole sector, which is the most relevant for Raman 
CDE processes. States with more than 5 % contribution to the energy-weighted sum rule 
(p!9|) are represented in the figure as triangles. They will be called "collective" excitations. 
They form three well defined bands accounting for, approximately, 7, 35 and 5% of the sum 
rule. The rest of the monopole strength is divided among 200 states with excitation energy 
lower than 30 meV. Fig. ^ shows also the lowest monopole SPE (triangles plus dotted line). 
A complex pattern of probability transfer between colliding levels, as the magnetic field is 
varied, is reflected in Fig. ^ in the form of abrupt variations of the number of collective 
levels. The situation is similar to the behavior of the dipole strength in the biexcitonEl. 

Dipole and quadrupole collective levels and the corresponding SPE in these sectors are 
shown in Fig. ^. One sees that dipole excitations rule, lower than monopole and 

quadrupole collective CDE, but the quadrupole SPE are lower at i? = 1 T, as mentioned 
above with regard to Fig. ^ 

We show in Fig. ^ a few results following from the pp-TDA calculations. In the upper 
figure, the lowest hflint for the intermediate state with A^ + 1 pairs, in which the added 
pair has le + lh = 0, is drawn. This magnitude can be taken as the renormalization of the 
absorption edge due to the background of A^ electron-hole pairs. Let us stress that there are 
two main effects contributing to this magnitude. The first is the blue shift induced by Fermi 
statistics, i. e. the added pair should occupy higher HF single-particle states. The second 
is the red shift caused by Coulomb (attractive) interactions. As can be seen in this figure, 
for A^ = 42 the net result is a blue shift of 12 - 16 meV. The apparent kinks are signals of 
ground-state rearrangements as the magnetic field is varied. 

The lowest part of Fig. ^ shows the dependence on A^ of the edge renormalization at 
i? = 1 T. It grows from 3 meV for 12 pairs up to 14 meV for the 42-pair system. That is, 
at a rate of 0.3 meV per pair in the dot. This magnitude can be used as a complementary 
way of determining the mean number of pairs in the dot. 

C. Raman spectra in CDE channels 

Let us consider Raman processes in which the final states are CDE. The first important 
question we would like to address is the role played by collective and SPE in resonant Raman 
amplitudes. 

We show in Fig. ^ the Raman differential cross section, computed from 
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= /^/^^it E m + -Ef- huj), (29) 

in which dAy is the sohd angle element in the direction of the dispersed light. We will use 
a smearing of the delta function as follows: 



with a phenomenological T f = 0.5 meV. 

The spectra in Fig. |^ are computed under conditions of normal incidence {(pi = (j)f = 
0, only monopole final states are excited) and parallel light polarization. The latter is 
supposed to disentangle collective CDE modes from SDE in electronic qdots under non- 
resonant scattering|ll. The monopole strengths are also included in the figure for comparison. 

The upper figure shows results at i? = 1 T. There are always Raman peaks associated 
to the more collective CDE states, although their magnitude rapidly vary with the incident 
laser frequency. The overall behavior, which is apparent in this figure, is that low-energy 
weakly collective or SPE are favored at "extreme resonance", i. e. when hui is near 1574 
meV in this situation, whereas 30 meV above the effective band gap high-energy collective 
or weakly collective CDE give the strongest peaks. 

Notice that the maximum peak intensities under extreme resonance are reached for laser 
frequencies a few meV above the renormalized band gap. Thus, as the laser frequency moves 
above 1574 meV, the amplitudes corresponding to weakly collective or SPE initially increase, 
but further experience a sudden drop. 

In the lower part of Fig. 0, the spectra at i? = 2 T are drawn. We notice variations 
in the peak distributions as compared to the B = 1 T results. Notice also that, in both 
spectra, the strongest CDE state is not seen as a distinct peak for hiVi = 1600 meV. 

The qualitative conclusions to be extracted from Fig. |^ are thus the following: (a) 
Comparable Raman intensities for strongly collective and for weakly collective states (even 
for SPE at intermediate excitation energies), (b) A complex pattern of variations of the 
Raman intensities as the frequency of the incident laser or the magnetic field is varied. 
Distinct peaks are seen only in certain intervals of these magnitudes, and (c) A richer 
structure of the Raman spectra as compared with the charged quantum dots. 

The second important point to discuss, from the qualitative point of view, is the excita- 
tion of high multipolarity modes at non-zero momentum transfer. We notice that the light 
wave vector is around 8 x 10^ nm~^ in the present situation. The maximum momentum 
transfer is thus 1.6 x 10^ nm~^ in backscattering geometry. We show in Fig. ^the quadrupole 
spectra at -B = 1 T in the parallel light polarization configuration and momentum trans- 
fer equal to 8 x 10'^ nm~^ (0, = 0/ = 7r/6). Of course, in an experimental curve all the 
multipolarities come together. We separate the quadrupole spectra to simplify the analysis. 

First, we notice that quadrupole Raman intensities are 1/10 of monopole ones. On 
qualitative grounds, one expects quadrupole intensities of order N^^^^^^{qi\\Dy, where D ^ 90 
nm is the system diameter, and Ngtates is the number of intermediate states participating 
in the process. {qi\\DY provides a factor 10~^, but the number of states contributing to 
quadrupole processes is roughly three times the states contributing to monopole processes 
(intermediate states with excess angular momentum 0, +1, and +2 in the A/ = +2 case. 
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for example). Thus N^tatesi^iW^)"^ ~ 10~^. Second, we observe an asymmetry between the 
Al = —2 and A/ = 2 spectra. Most of the A/ = —2 peaks correspond to SPE or weakly 
collective states. The Al = 2 peaks, two or three times more intense, are concentrated around 
collective states, which strengths are more uniformly distributed. The most collective CDE 
state with AZ = 2 shows up as a distinct peak only in a thin range of frequencies. Other 
multipoles show similar behavior. 

Thus, the conclusions coming from Fig. ^ are the following: (a) The intensity of CDE 
Raman peaks with multipolarity I are proportional to (|/| + and (b) Negative-Z 

peaks correspond mainly to very weakly collective or SPE states. The peak associated to 
the most collective CDE state is well defined practically at any hui. On the other hand, 
positive-/ peaks are stronger and show a dominance of collective states. 

A third interesting question to be addressed is related to the modes excited when the 
dispersed light polarization is orthogonal to the polarization of the incident light. The 
results for monopole states at -B = 1 T are presented in Fig. ^. Under extreme resonance, 
we observe peaks associated to SPE modes with excitation energies lower than 12 meV. In 
particular, the lowest SPE at 4 meV is clearly distinguished. Raman signals due to CDE 
states are strongly suppressed in these conditions. 30 meV above extreme resonance, the 
dominant peaks are located at higher excitation energies. They correspond to SPE or weakly 
collective states. 

IV. CONCLUSIONS 

We have computed the Raman amplitudes for the light scattered of a qdot which contains 
42 e-h pairs. In an attempt to identify the states giving rise to the strongest peaks, we 
compared the Raman intensities with the multipole strengths. The result is that both 
collective and SPE states play important roles in Raman spectra. Their relative weight in 
the spectra is seen to strongly depend on the external magnetic field, the polarization of 
the scattered light and the frequency of the incident light. Taken in a wider context, this 
conclusion suggests caution when making an assignment to an experimental Raman peak, 
and urges for theoretical calculations in parallel to the experiments. 

The explicit construction of the wave functions for the intermediate states, always in the 
framework of mean-field time-dependent approximations, allows us to consider extreme as 
well as non extreme resonance conditions. In the same way, the formalism allows for any 
wave momentum transfer or any kind of light polarization. 

We can not presently consider SDE final states. The reason has been explained briefly in 
the text: the SDEs are "2p-2h" states, which can not be modeled by the RPA approximation 
adopted in this work. It shall be said, however, that none of the papers or ^ accounted 
for valence-band mixing effects in the intermediate hole state in Raman SDE channels. 

However, on qualitative grounds, it can be argued that SDE final states shall give strong 
Raman peaks, may be even stronger than CDE states. The argument goes as follows. It 
may be seen that the factor determining the Raman amplitude is in fact the orbital factor, 
i. e. the overlapping between the electron and hole wave functions. In the symmetric, 
Ne = Nh, system we are studying, the overlapping is high in the intermediate states (both e 
and h above the Fermi levels), but low for CDE final states because one of the annihilated 
particles is above its Fermi level, and the other is below. For SDE states, however, both the 
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annihilated e and h are below their Fermi levels, and the overlapping may be high. Thus, 
Raman SDE amplitudes could be even stronger than CDE amplitudes. 

In the electronic qdots, Raman scattering in SDE channels goes through hole band 
mixing. Apart from the low overlapping in final states, one would expect the amplitude to 
be proportional to the light hole component of the hole wave function. Due to the fact that 
the Coulomb interaction is diagonal in the band indexes, a strong electronic background 
could depress valence band mixing, thus making SDE amplitudes even weaker. Research 
along this direction is in progress. 
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APPENDIX A: EVALUATION OF BAND AND ORBITAL FACTORS 

We give in this Appendix the expressions for the band and orbital factors entering Eqs. 
(§) and @. The ratio of the band factor Si ■ Pa-y to the magnitude iP, where P is the 
interband GaAs constant, is given in Tab. |. S'^ = ±1/2 is the spin projection over the z 
axis. 

Conventionally, we assign = —1/2 to the rrij = 3/2 electron state in the valence band. 
The e± components are defined as 

The band factor entering Eq. (|), i. e. e^* ■ p^^, can be also obtained from Table | if we 
replace by and viceversa. 

On the other hand, the orbital factor in Eq. (^) is computed from the HE one-particle 
functions, Eq. (|). Substituting (|^) into the expression for the band factor, Eq. @, and 
making use of the expansion: 

e''^'-"" ^ 1 + iqm ■ - ^{qm ■ f^f, (A2) 
where f\\ means the projection of r'onto the xy plane, we get 

J e^'?-^>L(f)0;,(f)dV ^ CTsd;: I {ks, is\i - ^do\h, -h) + 



s,t 

I ^{K,h\e-''^d, + e''^d_^\kt,-h) 



{ks, ls\e-^''^d2 + e^''^d.2\kt, -It)} . (A3) 

In this last equation, di are the one-particle multipole operators, which explicit expression 
is given in Appendix 0. The kg and Ig are, respectively, the radial and orbital quantum 
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numbers of the 2D oscillator state Xs- The angle Oi is by definition equal to zero, i. e. the 
z axis is oriented along 

The orbital factor entering Eq. (^) can be obtained formally from ( [A3| ) upon substituting 
i by / and taking the complex conjugate of the whole expression. 



APPENDIX B: MULTIPOLE MATRIX ELEMENTS AND SUM RULES 



The evaluation of one-particle elements, , requires the expansion (H) for HF functions: 



(Bl) 



s,t 



where the elements d[^, taken over oscillator functions, are given, when / = 0, by: 



{ks, h\%kt, It) = S{h, It) \ {2kt + \lt\ + mks, kt) - ^{kt + l){kt + \lt\ + l) 6{ks, kt + 1) 
to L 



- ^/kt{kt + \lt\) 6{ks,kt-l)^ 
Whereas, for Z > 0: 

, Min(;,fct) 

{ks,is\4\ktjt) = siisjt + i) E (-1) 



(B2) 



n 



ktl ikt-r + \lt + l\)\ 



iks-r + \lt\)\ 



r=o il-r)\r\\{kt + m {kt-r)\ 

x6{k„kt-r)] lt>0, 

Mm{l,ks) ^ 

x6{kt, ks - r); k < -I, 

Mm{l-\lt\,kt) 



r=0 



(l-\l,\-r)\r\{\lt\~ks + kt-r)\ 



X, 



ktl 



ks\ {kt-r + l)\ 



|2 



-e{ks- kt + r); 



-\\ ikt + \lt\)\ih + \lt + l\)\ ikt-r)\^ 

-I <lt< 0. (B3) 
Where 0(x) = 1 for < x < Min(|/f |, kg) and zero otherwise. Finally, for / < 0, we get: 

{k,,h\di\kt,lt) = {ktMd\i\\ks.lsY . (B4) 

On the other hand, the elements (A|(r^)^|A), entering the r. h. s. of the sum-rule 
equations (|l9| - pO| ) are evaluated as: 

{A{r'm = T.ClsCM{s\{ry\t), (B5) 



s,t 



where: 



2\« 



{ks,ls\ -p2 \kult) = S{ls,lt)^ 



ks\kt\ 



{k, + \is\m + \is\)\ h h ^ ^ 



x5{ks — m,kt — n) 



{kt-n+%\+i)\ 



— m)\m\{^ — n)\n\ {kt — n)\ 



(B6) 
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TABLES 



ce \ qh 
'-'za \ '-'2:7 


-1/2 


1/2 
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-1/2 








TABLE I. The quotient Si-pa^/{iP). 
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FIGURES 



FIG. 1. (a) Inelastic light scattering leading to final states which are CDE of the ground state, 
(b) An example of Raman scattering in SDE channels. 

FIG. 2. A set of electron and hole HF levels at B = 1 T. The Fermi energies are indicated as 
dotted lines. The less energetic transitions with AL^ = ±2 are represented by arrows. 

FIG. 3. Electron and hole densities in the HF ground state. 

FIG. 4. Monopole collective CDEs of the A/"-pair system and the lowest SPE. 

FIG. 5. Dipole and quadrupole collective CDEs and the lowest SPEs. 

FIG. 6. Absorption edge renormalization: (a) As a function of B for N = 42, and (b) As a 
function of A/" for B = 1 T. 

FIG. 7. Raman spectra under conditions of normal incidence and parallel polarization. Only 
final CDE states are considered. 

FIG. 8. Raman spectra of quadrupole CDE states, (pi = cpf = n/6 and parallel light polariza- 
tion geometry. 

FIG. 9. The same as in Fig. 7, but the polarization of the scattered light is orthogonal to the 
polarization of the incident light. 
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> 

>^ 

LU 



> 


C 
LU 




20-1 
15 
10 

5 


-5 
■10 
■15-1 
■20 



holes 



■ I ' I 

12 -10 -8 



I ' I ' I ' I ' I ' I ' I ' I ' 
-4 -2 2 4 6 8 10 12 



Angular Momentum 



CM 



0.009 

0.008 - 

0.007- 

E 0.006-1 
c 

0.005 - 

-I— > 

0.004- 

o 

Q 0.003 -I 
0.002 - 
0.001 - 
0.000 





0.009 - 




0.008 - 


Cvj ' 

1 


0.007 - 


E 


0.006 - 


H — ' 


0.005 - 


CO 







0.004 - 


Q 


0.003 - 






0.002 - 




0.001 - 




0.000 - 




r [nm] 



26 



J I I I I I I I I 



24- 



A 



22- 



A 



20- 



> 
CD 

E 18 



D) 

CD 
C 
LU 

O 
'■t—> 
C\3 



16- 



14 - 



o 12 

X 
LU 



A 
A 



A 



A 



A 



10- 



8- 



6- 



A SPE 
A Al = 



4- 



A- 



A 



A 



0,0 



0,5 



1,0 

B[T] 



1,5 



2,0 




> 

CD 

E 



18 



16- 



14- 



^ 12 



iS 10 

c 
o 

8-1 



o 

LU 



6- 



4- 



2 - 













► 






► 








► 






► 


► 


► 






< <j 






► 


Al = 2 






< 


Al = -2 






► 


SPE (+2) 






< 


SPE (-2) 


► 




-<i. 




















< 



0,0 0,5 



1,0 

B[T] 



1,5 



2,0 



> 

CD 

c 
o 

OS 
N 

OS 

E 

o 
c 

CD 

CD 
D) 
■D 
CD 

C 

o 



14- 
12- 
10- 



=1/2, S = 

e ' h 

O S =-1/2, S = 

e h 

I +1=0 

e h 

N = 42 



-1/2 

: 1/2 




B[T] 



1 



1 



1 



S = 1/2, S = 

e ' h 

I +1=0 

e h 

B = 1 T 



-1/2 



o 

CO 

< 



8 
6 
4 



10 15 20 25 30 35 40 45 

Number of pairs 




"E 



0.04- 



0.03- 



■t 0-02 

CD 



0.01 - 



0.00 



hv = 1580 meV 




Al = 


hv = 1600 meV 

V 




B = 2T 


Bars: Strength/1 W \J 


i J . I. 


.1 







1 

5 



10 15 

Raman Shift [meV] 



20 



0.004 - 



"E 

- I—' 

"(/) 
c 

CD 



0.003 - 



0.002 - 



0.001 - 



0.000 



Al = -2 




hv = 1573 meV 


B = 1 T 




hv = 1600 meV 

I 






Bars Strength/100 




,1 .... 


II.. , . h .. 



0.012- 



„ 0.010- 

^ 0.008 -1 

" 0.006 -I 
>^ 

S 0.004- 



hv = 1573 meV 

I 

hv = 1600 meV 

I 

Bars Strength/10 




10 15 

Raman Shift [meV] 



1 



Al = 2 
B = 1 T 



20 



Orthogonal polarization 



0.040 - 
0.035 - 
0.030 - 
0.025 - 
0.020 - 
0.015- 
0.010- 
0.005 - 
0.000 • 



hv = 1573 meV 

I 

hv = 1600 meV 



A! = 
B = 1 T 




"T" 
20 



10 15 

Raman shift [meV] 



